Correlation dynamics of strongly-correlated bosons in time-dependent optical lattices 
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We analyze by means of Matrix-Product-State simulations the correlation dynamics of strongly- 
correlated superfluid Bose gases in one-dimensional time-dependent optical lattices. We show that, 
as for the case of abrupt quenches, a quasi-adiabatic modulation of the lattice is characterized by 
a relatively long transient regime for which quasi-local single-particle correlation functions have 
already converged to a new equilibrium, whereas long-range correlations and particularly the quasi- 
condensate fraction may still present a very significant dynamics well after the end of the lattice 
modification. We also address the issue of adiabaticity by considering the fidelity between the 
time-evolved state and the ground-state of the final lattice. 
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I. INTRODUCTION 

Strongly-correlated atomic gases in optical lattices 
constitute one of the most active fields in the physics of 
cold gases. Nowadays spectacular progresses allow for an 
unprecedented degree of control of these systems, which 
permit the detailed analysis of many-body phenomena as 
e.g. the realization of the superfluid (SF) to Mott insu- 
lator (MI) transition in bosonic lattice gases [ij], the 3D 
fermionic Mott insulator 0, or the Tonks gas [1] . 

The high tunability and long characteristic time scales 
of these systems offer an ideal scenario to investigate non- 
equilibrium dynamics in a way not available in traditional 
condensed matter systems. In particular, lattice hopping 
rates may be easily tuned by modulating the intensity of 
the lasers creating the optical lattice, and the interactions 
may be also modified in real time by means of Feshbach 
resonances and time-dependent magnetic fields. These 
changes may be produced fast enough to be considered 
as a sudden quench. These quenches have attracted a 
growing attention in recent years, in particular in what 
concerns the evolution of correlations and possible equili- 
bration after a quench [1, S 0, i, H, [10, [li, 12] . Thcrmal- 
ization (or actually the absence of it) was recently studied 
in a milestone experiment performed in nearly integrable 
one-dimensional Bose gases by Kinoshita et al [l3| . 

On the contrary, if the modification of the system is 
very slow, much slowlier than the tunneling rate, one can 
in principle assume the evolution as adiabatic [14, 15, 1^. 
The issue of adiabaticity is however far from trivial, espe- 
cially in low dimensional gapless systems, as recently dis- 
cussed by Polkovnikov and Gritsev 17[ . Interestingly, in 



the so-called non-adiabatic scenarios, the adiabatic limit 
cannot be reached no matter how slow the change is in- 
troduced. Although this result is strictly speaking only 
applicable in integrable harmonic systems, it was shown 
in Ref. ^ITJ that this non-adiabatic scenario may be ob- 
tained by considering an initial non-interacting ID or 2D 
Bose gas in a lattice under a slow increase of the inter- 
action strength. The harmonic approximation remains 
accurate as long as Uq/Jtiq <^ 1, where (see below) Uq 
characterizes the on-site interactions, J is the hopping 



rate, and no is the filling factor. In this regime the trun- 
cated Wigner approximation (TWA), plus an additional 
first-order quantum correction, allows for an accurate de- 
scription of the correlation dynamics [13] • 



In this paper we study the correlation dynamics of a 
superfluid Bose gas in a ID lattice during and after a 
modification of the lattice depth. Contrary to the case 
of a quench 6, 7, 8, 9, 1 Q,, [U, [12] , the modification is 
not considered as instantaneous, but rather a finite linear 
ramp. In addition, contrary to the (quasi-)harmonic sce- 
nario discussed by Polkovnikov and Gritsev (l7j we are 
here particularly interested in the correlation dynamics in 
the deeply quantum regime C/q 3> J at low filling no < 1, 
where the system remains superfluid although strongly 
correlated. In this regime, quantum fluctuations are 
dominant, and hence TWA approximation cannot be em- 
ployed to describe the dynamics. To this aim, we employ 
time-dependent Matrix-Product-state techniques, which 
allow us to study accurately relatively large systems. We 
show that, as for the case of abrupt quenches the correla- 
tion dynamics is characterized by an intermediate regime, 
in which quasi- local correlations have converged to a new 
equilibrium, although long-ranged correlations, and in 
particular the quasi-condensate fraction still presents an 
observable dynamics (which keeps evolving well after the 
end of the slow ramp). Since the system is not integrable, 
eventually a new equilibrium is reached. We analyze this 
final state and the adiabaticity of the modification by 
means of the transferred energy and the fidelity with re- 
spect to the expected ground state solution. 



The scheme of the paper is as follows. Sec. |TT] intro- 
duces the model under consideration and the numerical 
methods employed. Sec. IIIII is devoted to the analysis 
of correlation functions and quasi-condensate fraction. 
Sec. IIVI studies the fidelity of the final evolved state with 
respect to the ground state of the final Hamiltonian. Fi- 
nally, Sec. fVl summarizes our conclusions. 
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II. MODEL AND METHODOLOGY 

In the following we consider bosons in a deep lattice 
constrained to the lowest energy band. In this regime, 
the free energy of the system is described by the BHH 
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where hi (a|) is the annihilation (creation) operator of 
a boson at the «-th site, hi is the corresponding num- 
ber operator, ^ is the chemical potential, J{t) is the 
time-dependent tunneling amplitude between neighbour- 
ing sites, and JJq > is the repulsive on-site interaction. 
We consider in the following that dX t — Q the system is 
in the ground-state for the initial Ji — J(0). 

As mentioned above, the system may be driven out of 
equilibrium either by modifying the hopping rate (as we 
consider here) or by modifying the interaction rate (e.g. 
by means of Feshbach resonances). In the following we 
consider a time-dependent lattice depth, which leads to 
a linear-ramp J(t) = Ji + {Jf — Ji)j- for an initial time 
interval < t < i^, where Jf is the final hopping. For 
t > tr the system evolves at a constant J = Jf. We 
consider sufficiently large post-ramp times such that the 
quantities of interest enter into a new equilibrium. 

We consider in our calculations a lattice with L = 60 
sites with open boundary conditions, which is sufficiently 
large to minimize finite-size effects at the lattice center, 
where we evaluate the correlations discussed below. In 
our time-evolution simulations, we work in the canonical 
ensemble with two different total number of particles, 
N = 20 which leads to an average lattice-site filling h ~ 
0.3 below half- filling (HF), and N = 50 which leads to 
n ~ 0.8, i.e. above HF. In all simulations discussed below 
we considered J, = 0.1375C/o and Jf = 0.2500C/o- These 
values were chosen relatively close to each other to allow 
for the convergence to a new equilibrium discussed below 
within a numerically available evolution time. In spite of 
that these values are sufficiently different to allow for the 
study of the ramping adiabaticity. Note that the hopping 
rates Ji and Jf are rather low and comparable to the 
critical tunneling for on-set of MI (tip of the lowest MI 
lobe), which for ID is found at J ~ 0.2Uo. In that regime 
quantum fluctuations are highly relevant, but due to the 
low filling h considered, the system remains within the 
(highly-correlated) SF regime. 

In our calculations we first obtained the ground-state 
at J — Ji by means of an MPS-algorithm using a similar 
approach as that of Refs. [l^, In this DMRG-like 
technique, the many-body state is approximated by a 
MPS ansatz of the form 
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where the A^s are matrices of a given dimension D and 
d is associated to the on-site dimension, which in our 



FIG. 1: Time-evolution of several correlation functions con- 
sidering a ramp-time of (a) t,- — 0.2 J/ and (b) tr — 10.0 J/. 
The used density population is n — 0.3. 



case represents the maximal number of atoms per site 
considered, (for our parameters it can be safely assumed 
as c? = 2). The MPS-algorithm performs in an efficient 
way a variational method using the matrix elements A^' 
as variational parameters, and basically reduces to sub- 
sequent DMRG sweeps over the lattice until achieving 
energy convergence [20| . 

Once the ground state is obtained with J ~ Ji we 
evolve the problem in real time for i > by means of 
the time-evolving block decimation (TEBD) method [2l|. 
TEBD takes advantage from the fact that the Hamilto- 
nian can be written as the sum over an even and an odd 
sites and then, the time evolution operator is approxi- 
mated using the Trotter-Suzuki expansion formula [23 |. 
Since the entanglement of the system grows on time, to 
keep the matrices of the same size means to loose infor- 
mation through out the time evolution. In order to avoid 
that, we adapt the value of D in each step to keep the 
state as faithful to the real one as possible. 

III. CORRELATION DYNAMICS AND 
QUASI-CONDENSATE FRACTION 

In the following we are particularly interested on 
the dynamics of the single particle correlation functions 
G{A;t) = (ajflj+A)- In our numerical calculations we 
evaluate G{A,t) at the lattice center j = 0. In equilib- 
rium, the ID SF phase is characterized by a polynomial 
decrease G{A) ~ A"'' with the relative distance A [2^ 
(contrary to the gapped MI phase where an exponential 
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FIG. 2: Spatial correlation for the initial state |(/5igs), the 
ground-state \^sga) of the final configuration and the evolved 
state \tpevoi) after t/J — 35 for a ramp time of tr — 10. OJ/. 



decay is expected). 

We have analyzed the evolution of G{A, t) / G( A, t — Q) 
for different values of A and different ramping times t^. 
Fig. [ija) exemplifies the case of an abrupt ramp (basi- 
cally an instantaneous quench) with tj. — 0.2 J/, whereas 
Fig.[TJb) depicts typical results observed for a mild ramp, 
in this case tr — 10.0 J/. Both cases are calculated at a 
filling factor n = 0.3. Note that G'(A,t) shows in both 
cases a significant dynamics following the linear ramp 
(due to number conservation the density, i.e. A = 0, is 
unaffected by the lattice modulation). 

Both abrupt and slow ramps lead to an evolution of the 
correlations characterized by an initial short-time scale, 
followed by an eventual convergence into a new equilib- 
rium at longer times (observe, however, that the cor- 
relation dynamics following the abrupt quench presents 
a short-time modulation which persists well within the 
quasi-equilibrium region ^). The short-time evolution 
of the correlations continues well after the end {t — tr) 
of the ramp (even for such a mild ramp). Notice also 
that short-distance correlations, in particular A = 1 con- 
verge significantly quicker than correlations at larger dis- 
tances. As a consequence the lattice bosons experience 
a transient regime characterized by a quasi-equilibrium 
of local or quasi-local observables coexisting with out-of- 
equilibrium global properties. After this transient regime 
the system reaches a final equilibrium, characterized by 
an equilibrium correlation G(A) (Fig. [2]). 

The transient regime becomes particularly clear from 
an analysis of the quasi-condensate fraction. This frac- 
tion may be defined as the largest eigenvalue Aq of the 
density matrix (ajflj) of the system '23|, and hence may 
be considered a global property of the system, influenced 
by correlations at any available A. Although strict con- 
densation is prevented in ID, quasi-condensation, charac- 
terized by a distinct finite Aq, is possible in finite systems. 

The time evolution of Ao(i)/Ao(0) is depicted in Fig.O 
Note that the quasi-condensate fraction evolves at a 
much longer time scale (much larger than the ramp time 
even in the mild-ramping case) than quasi-local correla- 
tions (it has not yet fully converged in our typical cal- 
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FIG. 3: The evolution of the quasi-condensate fraction for 
two different particle densities n = 0.3 and n = 0.8 in a mild 
ramp {tr — 10. OJ/). The horizontal lines denote where the 
quasi-equilibrium (A''^) values tend in each case. The dashed 
lines fit the initial curve stretch and the intersection points are 
defined as a qualitative estimation of the Aq variation time, 
ti{n = 0.3) = 17.9 and t2{n = 0.8) = 12.9. The plot is in 
units of Xoga where Xogs is the quasi-condensate fraction of 
the initial ground-state. 

culations). The filling factor n is of course important 
to determine the time scale of variation of the differ- 
ent correlation functions. The figure shows for a ramp 
tr — 10 J/ the evolution of the quasi-condensate fraction 
for the case of an average filling factor n = 0.8 com- 
pared to the case with n = 0.3. At higher densities the 
quasi-equilibrium is reached faster for the same ramping 
(a similar behaviour is observed for G{A)). As in Fig. [3] 
we may define a typical time scale for the variation of Ao, 
which is t{n = 0.3) = 17.9 and t{n = 0.8) = 12.9. 

IV. FIDELITY AND FINAL ENERGY 

A good tool for the analysis of the adiabaticity of the 
ramping is provided by the fidelity F = \{ip fgs\(p{t))\'^ , 
i.e. the Hilbert-space distance between the time-evolved 
state l'p(t)) and the expected ground-state \ipfgs) calcu- 
lated with the final J — Jf. The fidelity F is an inter- 
esting figure of merit for the adiabaticity of the ramping, 
since contrary to other figures (as the correlation func- 
tions discussed above) it just evolves while the ramping 
is on, since after the ramping the eigenstates of the final 
Hamiltonian are of course stationary, and consequently 
F remains constant. Hence, although other quantities 
require a rather long waiting time for comparing the fi- 
nal state and the time-evolved one, F provides an answer 
at the relatively short-time scale tr (see Fig. |4]). As ex- 
pected the abrupt ramping tr/ Jf — 0.2 leads basically 
to an instantaneous projection of the initial ground-state 
\^fis) into |<^/gs) (note that the overlapping is already 
rather large, 88%, due to the relative close values of Ji 
and Jf). As expected the milder ramping approaches 
further to Iv'/gs), however, the fidelity is still 5% off from 
\ipigg). Interestingly, this indicates that even very large 
ramping times significantly larger than the hopping time, 
and for a relatively small variation of A J = 0.1125?7oj do 
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FIG. 4: Fidelity of the time-evolved state to the ground state 
at the end of the lattice modulation for tr/Jf = 0.2, tr/Jf = 
10.0, tr/Jf — 20.0 and tr/Jf — 30.0 time modulation. 



not guarantee a perfect transfer into the ground state of 
the final configuration. The analysis of F for even milder 
ramps shows (see Fig. |4]) that milder ramps lead indeed 
to more adiabatic transfers (contrary to what may be 
expected in the harmonic regime p^j). 

Hence, although, as mentioned above, the correlations 
G(A) approach at longer times to a new equilibrium, 
this new equilibrium is not that given by the expected 
ground state, but by a new distribution with a higher en- 
ergy. The analysis of the final energy after the ramping 
does not provide however an equally strict adiabaticity 
analysis as that of the fidelity, especially for situations 
as those discussed here, in which Ji and J/ possess rel- 
atively close values. In our case the difference between 



the energy of \^Pigs) (or equivalently that of the system 
after the abrupt ramping) and that of is less than 

1%, whereas the time-evolution with the mild ramping 
provides a final energy less than 0.1% off that of 



V. CONCLUSIONS 

We have analyzed by means of MPS techniques the 
dynamics of correlation functions and quasi-condensate 
fraction of ultra cold lattice bosons in the deeply corre- 
lated superfluid regime during and after a finite linear 
ramp modulation of the hopping rate. We have shown 
that the evolution is characterized by a transient non- 
equilibrium state in which quasi-local correlation func- 
tions have already converged into a new equilibrium 
whereas long-range correlations and the quasi-condensate 
fraction present still a significant time dependence. Ad- 
ditionally, we have analyzed the formation (at a longer 
time scale) of a new equilibrium from an initial gas at zero 
temperature. By considering the fidelity with respect to 
the ground state of the final configuration we have shown 
that even rather mild ramps do not fully guarantee a per- 
fect loading of the new ground state. We have however 
shown that contrary to the harmonic regime p7| progres- 
sively milder ramps lead to a more adiabatic transfers. 
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